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1. 


INTRODUCTION : 


The  purpose  of  this  note  is  to  propose  calculation  procedures 
for  temperature  when  portions  of  a  lubricating  film  are  cavitated. 
These  procedures  are  suggested  for  use  in  conjunction  with  the 
treatment  of  temperature  as  presented  in  a  recent  paper  [Elrod, 
1991] .  They  are  to  be  coupled  with  a  cavitation  algorithm  [Elrod, 
1981] .  The  proposed  relations  implement  the  Jakobsson-Floberg- 
Olsson  cavitation  model  [See  Dowson  &  Taylor,  1979] .  That  model 
considers  the  cavitated  zone  to  be  one  of  constant  pressure, 
partially  occupied  by  striations  of  liquid  that  extend  rotationwise 
across  the  gap  in  Couette  flow.  Although  known  to  have  deficien¬ 
cies,  the  JFO  model  represents  a  considerable  improvement  over  ear¬ 
lier,  more  approximate  models.  Its  accuracy  will  be  tested  to  some 
extent  by  deviations  in  temperature  found  between  experiment  and 
the  present  analysis. 

Figure  1  shows  a  cross-section  of  the  central  portion  of  a 
submerged,  liquid- lubricated  journal  bearing.  This  particular 
configuration  was  chosen  because  it  involves  all  aspects  of  the 
proposed  analytical  modelling,  not  because  it  is  proposed  for  its 
practicality. 


Figure  2  Stretched  View  of  Lubri¬ 
cating  Film  Showing  Cavitated  and 
Full -Film  Regions 
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There  are  two  separate  problems  to  be  considered  here.  First, 
we  have  to  choose  an  appropriate  differential  equation  for  the 
temperature  within  the  cavitated  zone.  Second,  we  have  to  establish 
inter-relations  across  the  loci  of  separation  and  reformation. 

2.  the  temperature  EOTTATION  for  the  cavitated  Zp.NE: 


JFO  theory  presumes  that  the  liquid  in  the  cavitated  zone  is 
completely  in  the  form  of  striations  (e.g.  Figure  2) .  Accordingly, 
the  tPTTipP^rature  equation  for  this  zone  is  the  temperature  equation 
for  the  striations.  We  begin  by  considering  the  relations  for  a 
single  striation,  and  then  generalize  the  results  to  a  zone  defined 
by  separation  and  reformation  loci . 


Figure  3  Coordinates  Used  in  Analysis 

Figure  3  defines  the  coordinates  used  in  this  analysis.  The 
(iif f 0rent ial  equation  for  the  temperature  of  a  full  film  is  [Elrod, 
1990,  eq.  22]  : 

3T/5t  +  V*VT  +  {dT/dq)Dq/Dt  =  (4/h2)  O^T/9g^)  +  cp/pCp  [2.01] 

where : 

Dg/Dt  =  (-1/h)  [  (1+g)  (ah/3t)  +  V*hfjvdq]  [2.02] 

Here  D/Dt  is  the  Lagrangian  time  derivative,  the  vector  y  is 
formed  from  the  x  and  y  components  of  the  fluid  velocity,  h(x,t)  is 
the  film  thickness,  g  =  2z/h  -  1  and  cp  =  ii[  {du/dz)^  +  {dv/dz)^]  .  The 
fluid  density,  p,  and  the  specific  heat,  Cp,  are  taken  as  constant. 
The  operator  V  represents  ^(B/dx)  +  ^yid/By)  ,  with  both  spatial 

derivatives  taken  at  constant  g.  Thus  the  last  term  in  [2.02]  can 
equally  well  be  written  as : 
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becomes^ : 

o/ay)  [^ph^ap/ay) ]  =  i2(ah/at)  +  6u(a/ax)  [hd  -  ii/sio)]  [2.04] 

Let  US  assume  no  transverse  temperature  variation,  an  assumption  to 
be  be  verified  a  posteriori.  Then: 

(a^p/ay^)  =  12  (ah/at)  +  6u(ahj,/ax)  [2.05] 

Now  the  striation  is,  itself,  a  segment  of  full  film.  We  wish  to 
evaluate  Dg/Dt  as  determined  by  short -bearing  theory.  First  we  need 
dv/dy.  The  transverse  momentum  equation  gives: 

{d/d(;)\i{dv/dq)  =  (h/2)^ap/ay  [2.06] 

The  fluid  viscosity, p,  is  temperature  dependent.  Then  with  ^  =  l/p 
we  have: 

dv/dq  =  (h/2)  ^ (ap/ay) +  3^  [2.07] 

Or: 

V  =  (h/2)^  (ap/ay)  [2.08] 

At  the  upper  plate  the  transverse  velocity  vanishes.  Hence: 

0  =  (h/2)^(dp/dy)f^^gdq  +  [2.09] 

Therefore : 

v=  (h/2)^(ap/ay)  [j^^^^gdg  -  (lx/3lo)  f  ^dq]  [2.10] 

Take  the  y-derivative  of  eq.  [2.10],  neglecting,  as  before,  the 
temperature  variation  in  that  direction.  Now  substitute  for  d^p/dy^ 
its  expression  from  [2.05].  We  find: 

(a/ay)  (hv)  =  3[ah/at  +  (u/2)  (ahj,/ax)]^i„t/^p  [2.11] 

where : 

5int  =  [2.12] 

Note  that  if  there  is  no  transverse  temperature  variation  within 
the  striation,  then  there  is  no  transverse  variation  of  dv/dy. 

For  the  shear-driven  x-wise  velocity  component  we  have: 

u  =  (U/2Io)/^^^dg  [2.13] 


^  ^p  is  the  mean  fluidity  for  the  local  fluid  layer,  and  is 
given  by:  =  io  +  0.4^2  -  IiV3Io-  The  li  are  the  coefficients  of 
an  expansion  for  ^ (g)  in  Legendre  Polynomials. 
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[2.10] 


v=  (h/2)2(ap/ay)  [|^^^gdc  -  (i.i/3io)/_^^^dg] 

Take  the  y-derivative  of  eq.  [2.10],  neglecting,  as  before,  the 
temperature  variation  in  that  direction.  Now  substitute  for  3p/3y 
its  expression  from  [2.05] .  We  find: 

0/5y)  (hv)  =  3  [3h/at  +  (U/2 )  Oh^/ax)  ]  [2.11] 


where : 

-  i,/3i,/_'5<Jc  12-12] 

Note  that  if  there  is  no  transverse  temperature  variation  within 
the  striation,  then  there  is  no  transverse  variation  of  dv/dy. 


For  the  shear-driven  x-wise  velocity  component  we  have : 

u  =  (U/2io)  '• 

From  eqs.  [2.11]  and  [2.13]  we  then  obtain: 

V*hY  =  0/5X)  (hu)  +  0/ay)  (hv)  =  (U/2)  (3/ax)  [  (h/io) /^^^dg] 

+  3  [ah/at  +  (u/2 )  ah^/ax]  (^int/^p)  ^ 


where  h^  =  h(l  -  i.i/3i.o)  . 

To  determine  the  indefinite  integral  of  [2.14]  we  make  us  of 
the  formula : 

f‘^£^f(e)dedu  =  fUi-e)f(e)dG  [2.15] 

We  then  obtain: 

V.J_^^hydg  =  (U/2)  (a/ax)  (h/io)  [g/_^^^dg  -  f^g^dg]  + 

3  [ah/at  +(u/2)  (ahe/ax)  ]  (g-e)  (^i„t(6) /^p)de  [2. is] 

To  check  on  this  formulation,  we  carry  out  the  designated 
integrations  from  -1  to  +1.  The  result  is: 

J^V»hYdg  =  -2  (ah/at)  [2.17] 
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so  that  Dt^/Dt  firom  ag.  [2.02]  vanishes;  as  it  should,  at  both  the 
lower  and  upper  film  surfaces. 

In  the  case  of  constant  viscosity,  eq.  [2.16]  becomes: 

Dg/Dt  =  [(l-g2)/h]  [(g/2)  Oh/at)  +  UOh/ax)  (g+1) /4]  [2.18] 

Returning  now  to  the  differential  equation  [2.01]  for 
temperature,  we  note  that  if  5T/ay  =  0,  there  are  no  terms  in  the 
equation  to  generate  that  component.  Then  identically  the  same 
differential  equation  applies  across  the  striation  width,  and  the 
macrni tude  of  the  terms  is  independent  of  ths  striation  widbh. 
Although  the  above  analysis  was  constructed  for  diminishing  film 
thickness,  this  last  fact  encourages  us  to  use  the  relations  also 
for  a  divergincT  film  thickness .  The  logic  of  such  use  lies  in  the 
facts  that  (a)  films  are  known  to  support  some  tension  and  (b) 
rupture  into  a  multiplicity  of  striations  would  not  alter  the 
applicable  differential  equation.  We  propose  then  to  use  [2.02]  and 
[2.16]  to  trace  the  temperature  throughout  the  cavitated  zone,  from 
locus  of  separation  to  locus  of  reformation. 


3.  TR-RATMENT  OF  THE  INTERFACES : 

We  turn  now  to  treatment  of  the  interfaces  between  the 
cavitated  and  full-film  zones.  Figure  4  shows  a  small  section  of  a 
reformation  front.  The  flow  in  the  cavitated  zone  is  entirely  in 

the  x-direction.  Matching  flows  normal  to  the  front,  we  get: 

fficav  *  B  ~  ffiful  •  B  ~  [  P  '  ful  “  P  '  cav^  ^  ^  .01] 


Here  m  is  the  lineal  mass 
flow  per  unit  time,  p'  is 
the  film  mass  content  per 
unit  area  and  Vj,  is  the 
front  velocity  normal  to 
itself.  For  simplicity, 
we  now  write  this  rela¬ 
tion  for  a  constant -prop¬ 
erty  fluid  and  no  motion 
of  the  front. 


Figure  4.  Schematic  Diagram  of  Reformation 
Front 
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cos(a)©£^udz  =  cos  (a)  [Uh/2  -  (h^/l2/i)  0p/9x)  ] 

-  sin  (a)  (h^/12/x)  (9p/9y)  [3.02] 

Here  ©  is  the  fractional  liquid  content  on  the  cavitation  side  of 
the  front. 

Now  the  pressure  is  constant  along  the  front.  Therefore: 

cos  (a)  (dp/dy)  -  sin  (a)  (Sp/Ox)  =  0  [3.03] 

Substitution  of  [3.03]  in  [3.02]  gives: 

©U/2  =  U/2  -  (hVl2M)  (3p/9x)  [l+tan^  (a)  ]  [3.04] 

To  second  order  in  angle  a,  mass  continuity  is  expressed  by  the  x 
component  of  mass  flux. 

Figure  5  is  a  simplified  cross-section  of  the  streamlines  in 
the  neighborhood  of  the  reformation  front .  On  the  cavitated  side 
the  flow  is  shear-driven  until  close  to  the  front,  at  which 
location  the  striations  rapidly  spread  transversely.  On  the  full- 
film  side  there  may  be  a  region  of  reverse  flow  caused  by  the 
adverse  pressure  gradient.  Complex  though  these  phenomena  may  be, 
use  of  the  cavitation  algorithm  has  proved  to  give  a  satisfactory 
approximation.  However,  to 
treat  the  temperatures,  as 
contrasted  with  the  mass 
flows,  a  more  detailed 
analysis  must  be  attempted. 


In  Figure  5  we  show  the 
flow  from  the  cavitated 
region  passing  over  the 
recirculating  flow  on  the 
downstream  side  of  the  front . 
It  is  reasonable  to  suppose 
that  over  the  small  distances 
involved  that  the  temperature 
associated  with  a  fluid 
particle  will  not  change 
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much.  Thus  across  the  front  we  hypothesize  that: 

T  =  T  (Q)  [3 • 05] 

where : 

Q(g)  =  ©J'^udg  in  the  cavitated  zone  and  [3.06] 

Q(g)  =  J^udg  in  the  full-film  zone.  [3.07] 

Let  us  turn  now  to  the  recirculation  zone.  The  velocities  as 
shown  in  Fig.  5  are  not  precisely  correct,  for  Reynolds  Equation, 
as  an  asymptotic  representation  of  Stokes  Equation,  does  not  yield 
the  normal  velocities,  w,  that  would  correspond  to  the  bends  in  the 
depicted  streamlines.  Nevertheless,  on  the  basis  of  the  velocity 
distribution  shown,  we  invoke  again  Eq.  [3.05]  .  Progression  upwards 
from  the  wall  at  g=-l  yields  negative  values  of  Q.  The  maximum 
negative  value  of  Q  is  reached  where  u=0.  Beyond  this  value,  the 
Q's  repeat  until  they  finally  become  positive.  We  write: 

If  Q<0  then  T(Q,  u>0)  =  T(Q,u<0)  [3.08] 

The  point  where  the  x-wise  velocity,  u,  becomes  zero  is  found  from 
the  following  expression  for  velocity,  taken  from  Elrod,  1990. 
Thus : 

u=  [U- (2i.i/3)Bj  J_^^^dg/(2io)  +  [3.09] 

where  B^=  (h/2 ) ^  (3p/ax)  .  And  this  velocity  component  vanishes  when 
[|_^^^Cdg/J_^^^dg  -  i.i/(3i.o]  [(h/2)22ioOp/3x)  =  -U  [3.10] 

In  the  case  of  a  constant-viscosity  liquid,  this  last 
expression  reduces  to : 

[(h/2)M3p/ax)/MU]  (1-g)  =1  [3.11] 

In  order  for  this  equation  to  have  a  solution  it  is  necessary  for 

(h/2)M3p/ax)//xU  >  1/2  [3.12] 
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4. 


NTTMRRICAL  EXAMPLES  : 


The  foregoing  analysis  has  been  implemented  numerically.  As  an 
example  involving  all  aspects  considered,  we  treat  a  bearing  with 
the  following  geometry  and  lubricant  characteristics . 

Diameter  =  0.25m;  Length  .2m 

h„33„  =  .0001875m;  Eccentricity  =0.85 

Viscosity  =  .03  Ns/m^;  Thermal  Diffusivity  =  8E-8  mV s 
Volumetric  Specific  Heat  =  1.75E6  J/ (m^  °C) 

Surface  velocity  =  39.27  m/s 
All  surface  temperatures  =  50  °C 


bution  for  Bearing  Computations 


In  the  computations,  the  periphery  was  divided  into  30  even  incre 
ments .  Figure  6  shows  the  chosen  film  thickness  distribution. 
Figure  7  shows  the  fractional  film  content  obtained  with  the 
■yigcogity  held  constant  at  the  listed  value.  Under  the  circumstan 
ces,  the  film  content  and  velocity  distributions  are  independent  of 
the  liguid  temperature .  Figure  8  shows  the  velocity  distributions 
just  before  and  after  reformation.  As  can  be  seen  from  Fig.  8,  the 
index  13  is  for  the  position  immediately  before  reformation, 
whereas  14  and  15  lie  just  after.  Note  that  station  14  shows  some 
^0-y-0rge  flow  in  the  region  -1<Q'<.5  The  backward  flow  in  this  region 
must  return  in  a  layer  of  positive  velocity  located  adjacent  to 
g=0.5.  In  the  case  of  constant  liquid  viscosity,  these  velocity 
profiles,  as  well  as  the  fractional  mass  content,  are  independent 
of  the  number  of  internal  Lobatto  points,  so  long  as  that  number  is 
greater  than  or  equal  to  2 . 
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-1.0  -0.5  0.0  0.5  1.0 


LOCATION  IN  GAP 

Figure  8 .  Velocity  Distributions  in 
Neighborhood  of  Reformation 


-1.0  -0.5  0.0  0.5  1.0 


LOCATION  IN  GAP 


Figure  9 .  Tenperature  Profiles  in 
the  Neighborhood  of  Reformation 


Figure  9  shows  the  temperature  profiles  on  the  bearing 
centerplane  in  the  neighborhood  of  reformation  for  the  case  of  9 
internal  Lobatto  points.  Note  the  temperature  reversal  as  q 
proceeds  from  the  lower  wall.  This  reversal  is  in  line  with  the 
streamline  formation  schematically  shown  in  Fig.  5.  For  less 
numbers  of  Lobatto  points  this  behavior  is  represented  with  less 
and  less  fidelity.  For  three  Lobatto  points  it  is  shown  not  at  all. 


Nevertheless,  for  this  case  in  some  average  sense  the  features  are 
preserved,  for  at  station  15  (Figure  10)  there  is  little  difference 
in  the  resulting  temperature  distributions  from  3  and  9  points. 


-1.0  -0.5  0.0  0.5  1.0 

LOCATION  IN  GAP 


-1.0  -0.5  0.0  0.5  1.0 

LOCATION  IN  GAP 


Figure  10.  Comparison  of  Ten5)erature  Figure  11.  Shifting  Temperature  Profiles 
Profiles  at  Station  15  for  3  and  9  within  the  Striations 
Internal  Lobatto  Points 


In  Figure  11  we  observe  the  shifting  temperature  patterns 
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within  the  striations  of  the  cavitated  region.  For  steady-state 
operation  and  constant -viscosity  lubricant,  we  have  from  eq.  [2.18] 

DC/Dt  =  UOh/ax)  (1+0  (1-0) /4h 

Hence  if  the  film-thickness  is  increasing,  the  fluid  tends  to 
towards  the  upper  plate,  and  vice  versa.  In  the  absence  of  film 
convergence  or  divergence,  the  temperature  profile  is  symmetric 
between  the  walls,  peaking  at  the  middle. 


-1.0  -0.5  0.0  0.5  1.0 


LOCATION  IN  GAP 


Figure  12  Velocity  Distributions  with 
Variable  Viscosity 


LOCATIO.V  IN  GAP 


Figure  13 .  Temperature  Profiles  in 
Neighborhood  of  Reformation  with 
Variable  Viscosity 


Figures  12-15  show  some  results  obtained  with  the  fluid 
viscosity  varying  according  to  the  formula: 

fi  =  0.03  exp{-0.01448(T-50°C)  }  Ns/m^ 


The  fluid  viscosity  diminishes  with  temperature,  causing  less 
overall  temperature  rise  (See  Figs.  13,  14  and  15)  .  The  velocity 
gradient  midfilm  is  greatest  because  there  the  viscosity  is  least 
See  Fig.  12)  .  The  reformation  takes  place  slightly  earlier  (at 
index  13,  instead  of  14)  .  Again  downstream  there  is  little 
difference  between  the  results  for  3  and  9  internal  Lobatto  points 
(See  Fig.  14)  . 
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LOCATION  IN  GAP 


Figure  14  Temperature  Profiles  at 
Station  15  for  3  &  9  Internal  Lobatto 
Pts.,  Variable  Viscosity 


Figure  15  Shifting  Temperature 
Profiles  within  the  Striations 


5.  CONCLUSIONSj:, 

Analysis  and  algorithms  have  been  developed  for  the  treatment 
of  lubricating  films  involving  cavitation  and  temperature  effects. 
For  purposes  of  illustration,  the  techniques  have  been  applied  to 
a  simple,  externally- flooded  journal  bearing.  The  work  is 
considered  to  be  a  consistent  extension  of  Jakobsson-Floberg-Olsson 
cavitation  theory.  Deviations  between  the  present  analysis  and 
experiment  may  therefore  lead  to  a  better  assessment  of  the 
validity  of  JFO  theory^.  Furthermore,  the  results  presented  here 
can  be  used  as  test  cases  for  comparison  by  other  investigators 
with  other  programs. 
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2  For  example,  it  is  anticipated  that  the  temperatures  in 
the  cavitation  region  will  be  over-estimated  because  JFO  theory 
neglects  the  liquid  that  may  be  attached  to  one  surface  alone,  and 
therefore  not  subjected  to  shearing. 
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TJnMF.NflLATURE : 


specific  heat  at  constant  pressure 
film  thickness 

"convective"  film  thickness  =  h (l-i.i/3i.o) 
pressure 

defined  in  eq.  [3.06] 
time 

temperature 

X- component  of  velocity 
rotational  velocity  of  shaft 
y- component  of  velocity 
total  velocity  vector, 

surface  velocity  vector.  In  this  case,  gxU 
z- component  of  velocity 
rotational  direction  coordinate 
axial  direction  coordinate 

direction  normal  to  film,  with  value  of  0  at  half-way  point 

angle  between  the  downstream-pointing  normal  to  the  reforma¬ 
tion  front  and  the  unit  vector  in  the  x-direction 
integration  constant,  eqs.  [2.07]  -  [2.09] 

dimensionless  cross -film  location,  ranging  from  -1  to  1 
fractional  liquid  content  of  the  film  (p  dissipation 
function,  also  dummy  variable  in  eq.  [2.15]  . 

/i[(9u/9z)^  +  (9v/9z)^] 
liquid  viscosity 
liquid  fluidity,  l/fi 

coefficient  of  n*^*^  Legendre  polynomial  in  expansion  for  ^ 
mean  fluidity  of  film,  defined  in  footnote  of  section  2. 
liquid  density 
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